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Abstract 



> 

o 

We examine the effects of the iepton number vioiating R parity odd superpotentiai, W = XijkLiLjE'j:, 
C****" , on singie production of fermion (charginos and neutralinos) and scalar (sleptons and sneutrinos) super- 

partners at leptonic colliders for center of mass energies up to 500GeV — ITeV. The probability am- 
plitudes for all the five 2 — > 2 body processes: Z+ZJ — > xf^m, Xi v m (Xi^m), lmW ± , u m Z°(i' m Z ), ^7(^7), 
and the decays branching ratios for the produced superpartners are calculated at tree level. The 



rates for all five reactions are proportionnal to X^jj where J = 1, 2 for e e + and (1 fj, + collid- 



> 

ers, respectively. A semi-quantitative discussion is presented within a supergravity model assuming 
ly-^ . grand unification of gauge interactions and universal (flavor independent) soft supersymmetry break- 

' ing parameters mo (scalars), m 1 f 2 (gauginos) at the unification scale. The predictions obtained 

for the total and partial rates show that the single production reactions have a good potential of 
observability at the future e~e + and n~ fi + supercolliders. For values of the R parity violating cou- 
0^ ' pling constant of order 0.05, the x*' productions could probe all the relevant intervals for tan/3 

and mo and broad regions of the parameter space for the fi (Higgs mixing) and mi/2 parameters 
(\fj,\ < AQQGeV, TO1/2 < 240GeV), while the v and I productions could probe sneutrinos and sleptons 
masses up to the kinematical limits (mp < 500GeV, m ; ~ < 400GeV). Using the hypothesis of a single 



dominant R parity violating coupling constant, a Monte Carlo events simulation for the reactions, 



ZjZy — » xflm> x\,2 v m, Xi,2^m, is employed to deduce some characteristic dynamical distributions of 
the final states. 

PACS: 12.60.Jv 13.10. +q, i4.8 O.Ly, 13.85. Hd 
rN i T98/061 ^ep-ph/980750^ .... Phys. Rev. D: DU6463 

1 Introduction 

Should R parity turn out to be an approximate symmetry of the minimal supersymmetric standard model, 
the truly quantitative tests of such a possibility would have to be sought in high energy colliders physics, as 
was first emphasized in [Q, 0, || . The great majority of the existing theoretical studies for the LEP or the 
Tevatron accelerators physics have focused on signals associated with the LSP (lightest supersymmetric 
particle) decays and certain rare decays of the standard model particles (gauge 0, El [|, or Higgs 
bosons or top-quark 
for inos decays [jlO| ^ 

proposals were made recently to explain the so-called ALEPH anomalous four-jets events juj on the 
basis of R parity violating decays of neutralinos or charginos |ll| [lTj , squarks jl8|, |l^| , sleptons pQ| 

or sneutrinos ]2l| produced in pairs through the two-body processes, e + e~ — ► or e + e~ — ► //. 

(See ^] for recent updates and lists of references.) 

Apart from precursor studies devoted to the HERA collider J23|, [2^] , little consideration was given 
in the past to single production of supersymmetric particles in spite of the potential interest of a discovery 
of supersymmetry that might be accessible at lower incident energies. The reason, of course, is the lack of 

* Supported by the Laboratoire de la Direction des Sciences de la Matiere du Commissariat a I'Energie Atomique 



J). A few experimental searches have been attempted for Z° boson decays 
and also in more general settings |l2, 13 1. Proceeding one step further, interesting 
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information about the size of the R parity odd coupling constants other than the large number of indirects 
bounds deduced from low and intermediate energy phenomenology Q. Therefore, for obvious reasons, 
the existing single production studies have rather focused on resonant production of sneutrinos, charged 
sleptons 0, |, |12[ |7|, |§ |§ ||, || or squarks g, ||, EM The interpretation of the anomalous 



high Q events recently observed at HERA by the ZEUS |33|] and HI ]3J| Collaborations, in terms of 
squark resonant production, has also stimulated a renewed interest in R parity violation phenomenology 
& 

The collider physics tests of supersymmetric models without R parity entail an important change 
in focus with respect to the conventional tests: degraded missing energy, diluted signals, additional 
background from the minimal supersymmetric standard model interactions and uncertainties from the R 
parity violation coupling constants compounded with those from the superpartners mass spectra. Our 
purpose in this work is to discuss semi quantitatively the potential for a discovery and the tests of 
supersymmetry with 2 — > 2 body single superpartner production. Although several order of magnitudes 
in rates are lost with respect to the resonant single production, one can dispose here of a rich variety of 
phenomena with multilepton final states non diagonal in flavor. Besides, one may also test larger ranges 
of the sneutrino mass since this need not be restricted by the center of mass energy value. Encouraged 
by the recent developments on R parity violation and by the prospects of high precision measurements 
at supercolliders |S6| , we propose to study single production at leptonic (electron and muon) colliders 
for the set of five 2 — > 2 body reactions, Ijlj — > X*^, Ijlj — > X°Vm (x°v m ), Ijlj — * i^ l W ± , Ijlj — ► 
v m LZ° (u m LZ°), l^lj — > v m Ll {pmhl) [J = 1, 2], in a more systematic way than has been attempted so 
far. We limit ourselves to the lowest inos eigenstates. Let us note here that precursor indicative studies 
of the inos single production reactions were already presented in p?] , |38| and that recent discussions 
concerning the reactions, e ± 7 — > e^v and e ± 7 — > l^v, where the photon flux is radiated by one of the 



two beams, were presented in [39 . We shall restrict our study to the lepton number violating interactions 
LiLjE^. in association with the familiar gauge and Yukawa couplings of the minimal supersymmetric 
standard model. The final states consist then of multileptons with or without hadronic jets. 

This paper contains four sections. In section ^|, we present the main formalism for superpartners 
production cross sections and decay rates. In section ^, based on the supergravity approach to supersym- 
metry soft breaking parameters, we present numerical results for the total rates and the various branching 
ratios in wide regions of the parameter space. In section [|, we show results for final states distributions 
of the processes, — > X P> X ^ X°v, obtained by means of a Monte Carlo events simulation, using 
the SUSYGEN routine [^o). In section^, we state our conclusions. 



2 General Formalism 



Five 2^2 body single production reactions may be observed at leptonic colliders. We shall use the 
following short hand notation to denote the associated probability amplitudes: 



M(Xa + O 
M(x°a + V m ) 

M(l- L + W + ) 
M{v m + Z) 
M{v m + 7) 



M(l~ 
M{lj 
M(lj 
M(lj 
M(lj 



Xa + Vm), 
Vrn 4 



Z), 
7), 



(1) 



where J — 1, 2 is a flavor index for the initial state leptons (electrons and muons, respectively), the index 
a labels the charginos or neutralinos eigenvalues and the index m the sleptons or sneutrinos families. 
Our theoretical framework is the minimal supersymmetric standard model supplemented by the lepton 
number violating R parity odd superpotential, W = | X^jfc ^ijkLiLjE'f.. This yields the sfermion-fermion 
Yukawa interactions, 



1 3 

L = - Xijk[viLekR£jL + e.jLe-kRViL + &kR v iB. e jL — (i —>j)]+h.c. 



(2) 



[*^J,fc]=l 



where the sums labelled by indices, i,j,k, run over the three leptons and neutrinos families with the 
condition i j following from the antisymmetry property, Xijk = —\jik-- 
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Mass Intervals 


Decays 


Final State 


Amll 


m ; ~_ > m^- (1) 
m ; -_ < m^- (2) 


•X~ C^fc 






mc, > rn^- (3) 
mo < m%- (4) 


•X _ ~~ * ~~ * W^fc 


(B)i+i+Zrij 


7+/- p + e - 


mr, ma > m;- > m-^o (5) 
m^- > mj > m^o (6) 

m-£- > m£i > m^o (7) 


•X~ X%^p -> IpVpViljh 

•X~ -> "pip ->Vp l pX° 

* ^p^p^i^j^k 

•X~ -» ^ P -»• [pt'pX 

* IpVplSiljlk 


(c)i+i-ip?g 


t m l 'p c ° ' 
7+7- p ±/=F 


> m^- > m^o (8) 
m^- > m§ > m-^o (9) 


•XT -> X°1pQp, -* %q P Viljh 
•X~ -> 9p9p, ^q P qpX° 

^ QpQp^i^j^k 


(D)l+lp?I? + 2jet 


7+ e + e - 


m^- > m-^o + mw (10) 


•X" -> X°VF~ -> W-J/ii/fc 


(E)i+iti?w-g 


7+P+ e - 
'm 1 - 'm 



Table 1: TTie allowed chargino decays for different relative orderings of the superpartners masses. The column fields give 
the mass intervals, the decay schemes, the final states corresponding to the process, Ijlj — * X~Qn> with a single dominant 
coupling constant Ayj, and the leptonic components of the final states in the case of a single dominant coupling constant 
A m ii [m = 2,3]. The notation If stands for missing energy associated with neutrinos. 



2.1 Production Cross Sections 

Each of the processes in eq. (jl|) has a charge conjugate partner such that the transformation between 
pairs of conjugate amplitudes can be formally described by applying a CP transformation to the S- 
matrix. The relationship is most easily described at the level of the amplitudes squared obtained after 
summation over the spins. Because of the simple action of CP on the initial state, lj(k)lj(k'), it can 
be seen that the amplitudes for the pairs of charge conjugated processes are related by the substitutions, 
k <-> k' and «-> A*- fe . The tree level probability amplitudes are easily calculated by inspection of the 
Feynman diagrams given in Fig|l[ The formulas for the amplitudes are consigned in Appendix |a|. A few 
observations are in order at this point. First, the same configurations of lepton flavor indices, namely, 
A m jj with J = 1, m = 2, 3 for e~e + colliders and J = 2, m — 1, 3 for colliders, occur in all cases. 

Second, the amplitude for right chirality slepton l m B_ production has not been included in the above list 
of formulas for the reason that this is proportional to the coupling constants An m which vanishes by 
the antisymmetry property, Ay ft = —Xjik- Thirdly, all five processus can appear only in a single helicity 
configuration for the initial fermions (assumed massless), corresponding to identical helicities, namely, 
either l\l~r\ or (recall that physical helicity for antiparticle is opposite to chirality). Lastly, we 

observe that the relative signs between the s,t and u channels contributions are dictated by both the 
structure of the interaction Lagrangian and the signs of the Wick contractions for fermions. The results 
for the spin summed squared amplitudes are given by somewhat complicated formulas which we have 
assembled in Appendix Al We have checked that our formulas for x° and X productions agree with the 
results provided in J3?], 38| and EIJ . 

2.2 Decays 

In order to exhibit the possible physical final states, we need to consider the decays of the produced 
supersymmetric particles, taking into account both the minimal supersymmetric standard model inter- 
actions (denoted RPC or R parity conserving) and the R parity odd interactions (denoted RPV or R 
parity violating) . 
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Figure 1: Feynman diagrams for the processes, ijJT — > X~'m (a), l^lj — > X°^m W> 'j'j — * ^rn,L^ + ( c )' 'j'j — * 
t> m r,Z° (d) and l^lj — > 7 fej. The circled vertex correspond to the RPV interaction, with the coupling constant 

^mjjt an d th- e arrows denote flow of momentum. 
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Mass Intervals 


Decays 


Final State 




ma < (1) 


•V m — > Iklj 


(A)/ fe -/|^° 




mt, > m^+ (2) 






c c 'm'm 


ma, > m^+ (3) 
ma > m^+ > mj (4) 


* ^m^j^i^k 


(C) 


e +r 


ma > m^o (5) 


•V m -> i/ m X° -> VmViljh 




e + e _ , 


ma > m^+ > m^o (6) 
ma > m^+ > mj > m^o (7) 


* ^m^pVpVil j^k 
* ^mMp^pX, * ^m^p^p^i^j^k 


(E)i+i-i±i?qz° 


i+l~e + e~ 

7 + 7- e ±7T 


ma > m^+ > m^o (8) 
ma > m^+ > m^ > m^o (9) 


* ^Vn ^ X~^wi ~^ X QpQp 
* ^mQpQp^i^j^k 
• i/ m > X ^rn — ^ ImQpQp 
^ ImQpQpX. ^ ^mQpQpViljlk 


(F) l"£if + 2jet 


Z~e+e~ 

'm c 'm 


ma > n%x+ > m^a + mw (10) 


•Vm, — > X + im — > ^mf^ 

-> l m W + v i l j l k 


(G) l-l£l?$W+Z" 


Z _ e + e~ 



Table 2: The allowed sneutrino decays for different relative orderings of the superpartners masses. The column fields give 
the mass intervals, the decay schemes, the final states corresponding to the process, Ijlj — * Z°u mi with a single dominant 
coupling constant Ajj^ and the leptonic components of the final states in the case of a single dominant coupling constant 
A m ii [m = 2,3]. The notation If stands for missing energy associated with neutrinos. 



A number of hypotheses and approximations, which we list below, will be employed in the evaluation of 
partial rates. 

1) Supersymmetric particles decays are assumed to have narrow widths (compared to their masses) and 
are produced on-shell with negligible spin correlations between the production and decay stages. 
This allows us to apply the familiar phase space factorisation formula for the production cross 
sections. 

2) Spin correlations are neglected at all stages of the cascade decays such that the branching ratios 
in single or double cascades can be obtained by applying recursively the standard factorisation 
formula. 

3) Sleptons belonging to all three families and squarks belonging to the first two families are degenerate 
in mass. Therefore, for a given decay process as, for instance, x~ —> v p l v , either all three generations 
will be energetically allowed or forbidden. Furthermore, flavor off-diagonal channels such as, l\ — » 
h + Z ,.,. are closed. 

4) The lowest eigenstates of neutralinos Xa an d charginos xt (a = 1) are excited in the cascade chains. 

5) All superpartners decay inside the detector volume. In the presence of broken R parity, the condition 
for electric charge neutral LSPs to decay inside the detector yields comfortable lower bounds of order 
A>10- 7 [§f§. 

6) Either a single RPV coupling constant is dominant in both the production and decay stages, or a pair 
of RPV coupling constants are dominant, one in the production stage (A m jj) and the other in the 
decay stage (Xijk)- The latter case with two dominant RPV coupling constants may be of interest 
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Mass Intervals 


Decays 


Final State 


A m 11 


m i- < m x- C 1 ) 


•l~ -> hvi 


(A)/-^+ 


e~ 


> m x - (2) 


• l m -> V m f -> VmlkVjVi 




e~ 


m ; ~_ , m,y > m x - (3) 






e + e~l~ 


to ; ~_ > m^- > rnp (4) 


^ ^ m X ^ ^rn^i^i 
^ ^m^i^j^k 






> m^o (5) 


• l m — > / m X° — > ImViljh 


(D)Z-|±I?#W+ 




TOj_ > m^- > m x o (6) 


•i m -> X~v m -> v m _X°l p D p 




l p 6 6 , 

'p ° "TO 


m ; -_ > m^- > ma > m x o (7) 


*^rn ~* X ^r/i * V m lpVp 
^m^p^pX ^ l^m^p^p^i^j^k 






mj- > m x - > rn x o (8) 


•i m -> XT^m -> v m x°q p q p 

* ^mQpQp^i^j^k 


(F) Z± i?=#W+ + 2jei 


e + e~, 

P ±7T 


mj- > m x - > rriq > m x o (9) 


-> v m q P qpX° -> v m q p q p Viljl k 






to ; ~_ > m x - > m x o + mw (10) 


•lm -> X~Vm -> ^m_X°W^" 

— > v m W~Viljl k 


(G) 


e + e~, 

c 'm 



Table 3: TVie allowed slepton decays for different relative orderings of the superpartners masses. The column fields give 
the mass intervals, the decay schemes, the final states corresponding to the process, Ijlj — * W + l^, with a single dominant 
coupling constant Ayj, and the leptonic components of the final states in the case of a single dominant coupling constant 
Amil [m = 2,3]. The notation If stands for missing energy associated with neutrinos. 



since strong bounds on quadratic products exist only for a few family configurations. The strongest 
bounds arise from the fi — ► 3e decay Q: X p iiX p i2 < 6.5 10~ 7 , A p2 iA p ii < 6.5 10 -7 [p = 2, 3], while 
other quadratic product bounds are of order 10 -3 , 10 -4 . Besides, as long as the coupling constant 
A, which controls the RPV decays, is small in comparison with the gauge coupling constants but 
not very much smaller (so that the LSP decays inside the detector), then the branching ratio will 
depend weakly on A since the last stage of the decay chain (LSP decay) is independent of A. 

7) The widths for the decays with four and higher body final states are neglected, such as those 
which occur in slepton (sneutrino) decays for, m x -,m x o > toj (rn x -,m x o > mp), mediated by 
virtual charginos or neutralinos, namely, l m — * v m lkVjVi and l~ m — > l-mhlj^i (y m — > l m hljh and 

8) A supergravity model for the soft supersymmetry breaking parameters is used where, generically, 
X? is the LSP. 

The consideration of the various order relations in the superpartners mass spectrum leads to a list 
of decay schemes for the initially produced superparticles. These are displayed in Tables |], || and ||, for 
X , vl and II, respectively. The signals for the Xi decays are very few in number and will be discussed 



separately in section 4.2.1. Some comments on these tables are in order. Except for hadronic dijet pairs 
from the decay processes, x — > x ??', all other final particles will consists of multileptons and missing 
energy associated with neutrinos. In the hypothesis of a single dominant RPV coupling constant in the 
decays stage, one can deduce the various final states flavor configurations by an inspection of the tables 
(cf. tables captions). The produced x° > X ± I P 1 will decay according to cascade schemes dictated 
by the superpartners mass spectrum. It is important to distinguish the direct RPV induced decays: 
XT -* ViVjlk, xT -> kljh, X° -> Viljh, X° -» vjjlk, Pi -> ijj, l~ R -> IjUi, and IJ L -> l k Vi, from the 
indirect RPC induced decays: 1^ — > x~v, vl — ► X + '> R ~ * X°h &l — * 5d v i X° ~* l~h X° — * 
X~ — > l~v, x~ ~^ lv and x* 1 ~^ x°^ ± f° r two-body final states and, x~ ^ X°ff (f=leptons or quarks), 
for three-body final states. All the formulas needed to evaluate the partial decay widths are quoted in 
the Appendix [b| As can be seen from Tables [l]|,|, a given final state can arise from different processes, 



G 



depending on the relative orderings of the masses. A reaction chain occuring through an intermediate 
particle which is produced on-shell leads obviously to the same final state when the production of this 
particle is kincmatically forbidden and it must then occur through a virtual intermediate state. In the 
approximation of family degenerate sleptons, sneutrinos and squarks, the index p in the tables runs over 
the three generations. A single exception is the hadronic decay, x — > x° u pdp(x d p Up), which is restricted 
to the first two families because of the large top-quark mass. 

Another subtle point concerns the multiplicity of a given signal, namely, the number of different 
configurations which can lead to the same final state. Due to the antisymmetry property of Ay^, the 
final states from chargino RPV decays (cf. A and B in Table |l|) have a multiplicity of two. The reason 
is that these decays proceed through the exchanges of the sleptons (or sneutrinos) in families i and j, 
for a given Ay/j. This fact is already accounted for in the virtual x three-body decays (A(1),B(1), Table 
|l|), but must be put by hand in the cascade decays proceeding through the on-shell production of 
sleptons or sneutrinos which decay subsequentially (A(2),B(2), Table Q). 

To get a better understanding of the interplay between RPC and RPV decays, it is helpful to note that 
the branching ratios can be written formally as, Bo = c ■> f° r direct decays and, B\ = — s g i^_yi ' ^ or 
two-stages indirect decays, where B = - g 2 +A 2 is the LSP branching ratio, A and g are symbolic notations 
for the RPV and RPC coupling constants and c, Co are calculable constants. Of course, B = 1, if the last 
decaying particle is the LSP, which is the generical case. For values of the RPV coupling constants small 
with respect to the gauge coupling constants (A < 0.05), namely, g 2 ^> A 2 , the dependence on A of the 
indirect decays branching ratios is weak and we have, Bi 3> Bd- For large enough RPV coupling constant 
(for example, A = 0.1) or for suppressed indirect decays (due for example to kinematical reasons), the 
direct decays may become competitive and both the direct and indirect branching ratios depend strongly 
on A. 



3 The model and its parameter space 

We shall develop the study of single superpartner production within a non minimal supergravity frame- 
work, assuming the existence of a grand unified gauge theory and of family universal boundary conditions 
on the supersymmetry breaking parameters. The renormalization group improved classical spectrum of 
the scalar superpartners is determined in principle by the full set of soft supersymmetry breaking param- 
eters at the unification scale, Mx, namely, tuq (common scalars mass), (common gauginos mass), A 
(trilinear Yukawa interactions), B^ (bilinear Higgs interaction); by the parameters tan/3 = ^ = <g"> 
and fi(t), where t denotes the running scale; and by the gauge coupling constants, g a (t), along with 
fermions masses, m 2 (t). If one neglects the Yukawa interactions of quarks and leptons with the Higgs 
bosons, then the running masses of all sfermions remain family degenerate down to the electroweak 
breaking scale where they are described by the familiar additive formula, 



mjit) = m}(t) + ml + c f {t)m\ /2 ± m| cos(2/3)(T/ - Q f x w ), (3) 

where c/(i) are calculable coefficients depending on the gauge interactions parameters and the last 
term represents the D-term contribution, the upper and lower sign being for the left and right chirality 
sfermions, respectively. The most relevant Yukawa coupling constants, namely, those of the third family 
of up-quarks or, for large tan/3, of d-quarks and leptons, are expected to induce downwards shifts for the 
third family squarks (up and down) and sleptons, which depend non trivially on the parameters A and 
/i. In this work, we shall restrict consideration to the simple case of family independent running masses 
and employ the approximate representation in eq. (0) with the numerical values quoted in fl44|. Note that 



the total rates do not depend on the squarks masses and, as already remarked in section 2.2, the third 
families of squarks are not considered in the cascade decays. The charginos and neutralinos classical mass 
spectra are determined by the subset of parameters: M\{t), Mz{t), fi(t) and tan/3. For fixed m 1 / 2 , the 
solution of the one loop renormalization group equations is given explicitly by, TO1/2 = (1 — f3 a t)M a {t), 

where t — log(^jg-), Q denoting the running scale, (3 a = Ij^f, b a = (3,-1,-11) with a = (3,2,1), 
corresponding to the beta functions parameters for the gauge group factors, SU(3), SU(2)l, U(1)y, and 
gx is the coupling constant at unification scale. Note that the wino and bino masses are related as, 
Mi(t) = |M 2 (t) tan 2 9w- It is useful here to comment on the relation of our framework with the so-called 
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A=5, tan (3 = 2 A = 5, tan (3 =50 



Figure 2: The solution n(tz), at scale mz, for the electroweak symmetry radiative breaking equations, at running top- 
quark mass, mt 

(m t ) = 171 GeV ( m P°' e = 180 GeV), is plotted as a function of 10 ^ eV and 100 q^ v for four values of the 
pair of parameters, A and tan/3. 
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Figure 3: Mass spectrum for the chargino X\ an d the first two lowest mass neutralinos, \i an d x'j, as a function of /u. 
Four choices of the parameters, tan fi and M2 ( in GeV ), are used, as indicated on top of each window. 
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minimal supergravity framework in which one assumes a constrained parameter space compatible with 
electroweak symmetry breaking. Let us follow here the so-called ambidextrous minimal supergravity 
approach p5[ , where one chooses [mo, m^, A, sign(fi), tan/3] as the free parameters set and derives 
fi(tz), B^(tz), at the electroweak symmetry breaking scale, t z = ln.M^/m%, through the minimisation 
equations for the Higgs bosons potential. For fixed mo, muz and tan/3, varying A will let the parameter 
/J,(tz) span finite intervals of relatively restricted sizes. In figure ^, we give results of a numerical resolution 
of the renormalization group equations which show the variation of |/i(tz)| as a function of m and m^, 
and also exhibit its dependence on A. Note that the equations admit the symmetry, fJ>(t z ) —* —\x(tz)- 
Observing that /J>(t z ) is typically a monotonous increasing function of A, we see from Fig.|| that the 
corresponding incremental increase, 5fj.(tz)/n(tz), as one spans the wide interval, A € [—5, +5], is small 
and of order 20%. 

In the infrared fixed point approach for the top-quark Yukawa coupling, tan 8 is fixed (up to the 
ambiguity associated with large or low tan/3 solutions) in terms of the top-quark mass, nit = C sin/3, 
with, C ~ 190 - 210GeV, for, a 3 (m z o) = 0.11 - 0.13 g§. The dependence on A of the electroweak 
constraint also becomes very weak, so that fi(tz) is a known function of mo, mx/2 and tan/3 ]46| : 



9 ul 7 ,l + 0.5tan 2 /3 9 0.5 + 3.5 tan 2 8 

a H — -f - = mi 5— - + mf /2 5— -. (4) 

p 2 tan 2 8-1 1/2 tan 2 8-1 ' 



In section 4.2, we will discuss results for the branching ratios in this constrained model. The total 
rates are not affected in any significant way by which version of the supergravity models is used, since, 
as we will see, their dependence on tan/3 and mo turns out to be smooth. 

The main uncertain inputs are the superpartners mass spectrum and the coupling constants \jk- To 
survey the characteristic properties of single production over a broad region of parameter space, we found 
it convenient to consider a continuous interval of variation for /i(tz), namely, n(tz) £ [— 400, +400]GeV r , 
while choosing suitable discrete values for the other parameters: M 2 (tz) = 50,80,100,150,200 GeV, 
m = 20, 50, 150 GeV and tan/3 = 2, 50. We shall set the unification scale at M x = 2 l0 16 GeV and the 
running scale at Q 2 — m 2 z . For definiteness, we choose the coupling constant, which controls the size 
of the production cross section, at the reference value: X m jj = 0.05. This is the strongest bound for a 
slepton mass of lOOGeV |^6[ . The dependence of integrated total rates on \ m jj is then given by a simple 
rescaling ( 005 ) 2 but that of branching ratios on Ay^ (which may or may not be identified with X m jj) 
is more complicated because of the interplay between the RPC and RPV contributions which add up in 
the total decay widths. The reference value used here, Xijk — 0.05, is also an interesting borderline value 
since below this value the dependence of branching fractions on Ay& becomes negligible in generic cases. 

It will prove helpful in the following discussion to keep within sight the spectrum for the low-lying 
inos. We display in Figure || the results obtained by solving numerically the eigenvalues problem for the 
charginos and neutralinos mass matrices. Recall the current experimental bounds p7| , m^o > 23GeV, 
m^± > A5GeV, m^ > 37.1GeV and mj > A5GeV. The following remarks about Figure [| will prove 

relevant for the discussion on branching fractions: (i) The symmetry of the spectra under, [i —ji, is 
spoilt at low tan/3 as can be seen on the explicit expression for the inos masses in ]48[| ; (ii) The mass 
differences x + — \° increase with with a steep rise appearing at, ji = M2, the borderline between the 
Higgsino and gaugino dominant regimes; (iii) The spacings X2 ~ X? an d Xi — X? decrease in magnitudes, 
relatively to the Xi mass, with increasing M%. Although we show here the results for \\ mass, the 
interesting possibility of exciting the second neutralino \\ is not considered in the subsequent discussion. 



4 Results and discussion 
4.1 Total production rates 

The total production rates are evaluated by taking the angular integral, a = f^ Xrn ^dx, [x = cos#], 
over the differential cross sections which are given explicitly in eqs.( A.2)-( [A.6| ) in Appendix |a|. To follow 



the usual practice we shall set an angular cut-off to account for the poor detection condition along the 
beam pipe: 170° > 9 m > 10°, corresponding to Xjfi — COS Urn. — 

0.9848. 
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Figure 4: The integrated cross sections for the process, l^lj — > X\^rn, at a, center of mass energy of 200GeV , are shown as 
a function of fj, for discrete choices of the remaining parameters: (a) tan f3 = 2, mo = 50GeV, (b) tan/3 = 50, mo = 50GeV, 
(c) tan/3 = 2, mo = 150GeV and fa!,) tan/3 = 50, mo = 150GeV, with X m jj = 0.05. The windows conventions are such 
that tan /3 = 2, 50 horizontally and mo = 50, 150Ge V vertically. The different curves refer to the values of Mi of 50GeV 
(continuous line), lOOGeV (dot-dashed line), 150GeV (dashed line) and 200GeV (dotted line), as indicated at the bottom 
of the figure. 
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Figure 5: The integrated cross sections for the process, I'jlj — > Xilm, a * a center of mass energy of 500GeV , are shown as 
a function of fj, for discrete choices of the remaining parameters: (a) tan /3 = 2, mo = 50GeV, ^ tan/3 = 50, mo = 50GeV , 
(c) tan/3 = 2, mo = 150GeV and (d) tan/3 = 50, mo = 150GeV, with X m jj = 0.05. The windows conventions are such 
that tan /3 = 2, 50 horizontally and mo = 50, 150Ge V vertically. The different curves refer to the values of Mi of 50GeV 
(continuous line), lOOGeV (dot-dashed line), 150GeV (dashed line) and 200GeV (dotted line), as indicated at the bottom 
of the figure. 
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Figure 6: The integrated cross sections for the process, Ijlj — > Xi&m, at a center of mass energy of200GeV, are shown as 
a function of fj, for discrete choices of the remaining parameters: (a) taxi fi = 2, mo=50GeV, (b) tan/3 = 50, mo = 50GeV , 
(c) tan/3 = 2, mo = 150Gey and (d) tan/3 = 50, mo = 150GeV, with X m jj = 0.05. The windows conventions are such 
that tan /3 = 2, 50 horizontally and mo = 50, 150GeV vertically. The different curves refer to the values of M2 of 50GeV 
(continuous line), lOOGeV (dot-dashed line), 150GeV (dashed line) and 200GeV (dotted line), as indicated at the bottom 
of the figure. 



13 



(a) 




-400 -200 200 400 



V EGgV] 
(c) 

9 ^ — . — , — . — , — ^- 




-400 -200 200 400 



V [GqV] 



(b) 




-400 -200 200 400 



V [GqV] 
(d) 




-400 -200 200 400 



V [GqV] 



50GgV — 150GgV 

— -IOOGgV — - 200GqV 



Fi gure 7'. The integrated cross sections for the process, l~jlj — y y^'ra, a ^ a center of mass energy ofSOOGeV, are shown as 
a function of fj, for discrete choices of the remaining parameters: (a) tan (3 = 2, mo = 50GeV, ^ tan/3 = 50, mo = 50GeV, 
(c) tan/3 = 2, mo = 150GeV and ('dj tan/3 = 50, mo = 150GeV, with X m jj = 0.05. The windows conventions are such 
that tan /3 = 2, 50 horizontally and mo = 50, 150GeV vertically. The different curves refer to the values of M2 of 50GeV 
(continuous line), lOOGeV (dot-dashed line), 150GeV (dashed line) and 200GeV (dotted line), as indicated at the bottom 
of the figure. 
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Figure 8: The cross sections for the processes, l^jlj — » lmW + (a), l^jlj — > v m Z° (b) and l^lj — > v m 7 (c), are shown as 
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considered are 200, 500 and lOOOGeV, as quoted in the top window. 



15 



4.1.1 Inos production 



The results for the integrated rates of the production of the lowest mass eigenstates \i an d Xii & t LEPII 
energies, are displayed in Figures || and ||, respectively. The inos production rates depend smoothly on 
tan/3, and on the mass parameters, /i, mo, M 2 , in a way which closely reflects on the mass spectrum. 
Thus, the symmetry under /i «-> —fx is upset only for low tan/3 and the rates decrease with increasing 
M 2 . The only cases where fast variations of rates arise are for values of m and M 2 at which the center 
of mass energy hits on the sneutrino s-channel pole, yfs = mo- As mo increases, the resonance occurs at 
smaller values of M 2 since the sneutrino mass depends on M 2 , mo and tan/3 (see eq.(|J)). The pole cross 
sections themselves, as parametrized by the conventional formula, 



8tts T{v -> l+l-)T(i> -> X) 
™% (s-m?) 2 +T? 

4 lO s ( 100GeV ) 2 B(t> l+r)B{v X) /barns, (5) 



can grow to values several order of magnitudes higher. For clarity, we have refrained from drawing 
the cross sections close to the resonant energy in the same plot. This is the reason why the curves 
corresponding to M 2 = IbOGeV do not appear in Figures [|(c) (d) and |(c) (d). The effect of the pole 
can be seen for M 2 — 200GeV in Figures [|(a)(b) and ^(a)(b). We note also that for fj, = 0, Xi is 
a pure higgsino and the \i production cross section vanishes. The results for inos production rates 
at NLC or /i + u~ colliders center of mass energies are displayed in Figures || and [?]. The drop with 
respect to the LEPII energies is nearly by one order of magnitude. The second neutralino production 
rates, cr(x 2 ) = a(ljlj — > X2 u m), when this is kinematically allowed, turns out to be of the same order of 
magnitude as er(xi)- For y/s = 500GeV, <j{x\) and <r(x 2 ) are numerically close throughout the parameter 
space of our model. However, for yfs = 200GeV, there are regions (large tan/3, u < 0) where one has 
"■(xi) ~ 2er(xi) and other regions (low tan/3, a > 0) where one rather has cr(x 2 ) ~ ^(Xi)- As for the 
production rate of the second chargino, cr(%2~), this is always nearly an order of magnitude below <j(xi)- 



4.1.2 Sleptons production 

The slepton and sneutrino production rates depend solely on the sleptons masses and X m jj. The results, 
obtained by setting mr — mo, are displayed in Figure || for three values of the center of mass energies. 
An account of the mass difference between rrij and mo would not change the numerical results in any 
significant way. 

The differential cross section for the reaction Ijll — > £7 must be treated with special care because 
of its extreme sensitivity at the end points, x = ±1, in the limit of vanishing electron mass, m e — > 0. 
As appears clearly on the expression of the squared momentum transfer variable, t = (k 1 — p') 2 = 
m 7 — " m I + m 7)(l ~ "^T'e" 2 ')' f° r m 7 = ^ ne t-channel amplitude has a collinear singularity, t — > 
as x — > 1. An analogous collinear singularity occurs for the u-channel amplitude, u = (k — p 1 ) 2 — » 
as x — > —1. Imposing the cut-off on the center of mass angle 9 makes the regularisation of collinear 
singularities pointless. 

In the limit of vanishing m 7 , independently of x and m e , the sneutrino production cross section 
becomes infinite at the limiting energy point, yfs = mp. This accounts for the property of the numerical 
results for the integrated cross section to rise with mp, as seen in Figure ||(c). However, if one were to 
set m 7 at, say, the p-meson mass, in line with the vector meson dominance hypothesis, one would rather 
find the opposite behaviour with respect to the dependence on mo- Observe that the increase of the 
cross section with mo corresponds to the fact that, for mo ~ \/s, the process Ijlj — * £7 behaves like a 
sneutrino resonant production, accompagnied by the initial state radiation of a soft photon. 



4.1.3 Discussion 



In summary, the single production rates range from several 10 's of fbarns to a few 100's of fbarns at 
LEP energies and several units to a few 10's of fbarns at NLC energies. Therefore, the superpartners 
single production are at the limit of observability for LEPII assuming an integrated luminosity per year 
of 200pb~ 1 at yfs = 20QGeV. The prospects for single production should be rather good at NLC [ff9 
and p + /i~ colliders (50) since the assumed integrated luminosity per year is expected to be about 
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at yfs = 500Gey. Moreover, it is important to note here that had we considered for the RPV coupling 
constants, constant values for the product X m jj{ ioogbV )' ra ther than for X m jj, the rates would get 
an important amplification factor ( iooGeV ) 2 ^ or increasing superpartners masses. Note that the slcpton 
involved in the bound is of right chirality and thus is of opposite chirality than the slepton involved in the 
rate. Of course, the masses of II and Ir are related in a given model. At y/s — 500GeV and assuming 
Kjk > 0.05, all the single production reactions should be potentially observable over a broad region 
of parameter space. The slepton production reactions could then probe slepton masses up to AOQGeV 
(Figure ||(a)) and sneutrino masses up to 500GeV^ (Figure ^(c)). The ino production reactions could 
probe a large region of the parameter plane (ji, M2), since the dependence on the parameters Too and 
tan (3 is smooth. To strengthen our conclusions, it is necessary to examine the signatures associated with 
the final states, which is the subject of the next section. 



4.2 Branching Ratios 

In the narrow resonance approximation, the partial transition rates are readily obtained by multiplying 
the total rates for each reaction with the decay branching fractions. The various final states for each 
of the 2^2 single production reactions have been listed in Tables [l], |^ and 0. The leptons family 
configurations in the final states will depend on the hypothesis for the RPV coupling constant (single or 
pair dominance). 

With the purpose of testing characteristic points of the parameter space, we have evaluated the 
branching ratios for the decays of the superpartners, namely, Xi and for variable /i at discrete 
choices of M2, too an d tan/3, such that the main typical cases in the ordering of the masses to^o, to~±, 
to ; ~, Trip , can be explored. The results are shown in Figures 0, [h] and O, for the chargino, the sneutrino 
and the slepton decays, respectively. The curves for the various branching ratios are distinguished by 
the same letters (numbers) as those used in Tables [j], || and || to label the various final states (decay 
processes). We shall now discuss in turn the various superpartner decay schemes corresponding to the 
five single production reactions. 



4.2.1 Lowest mass Neutralino 

The branching ratios for the \i desintegrations are best analysed separately. For convenience, we do not 
treat the cases, to^o > rriq and to^o > to ~± , since these arise marginally in most of the currently favored 
models (supergravity or gauge mediated soft supersymmetry breaking) . The cascade decays which occur 
if to~± < to^o are also not considered since the corresponding region of the parameter space (Figure |^(c)) 

is forbidden by the experimental constraints on the inos masses. Thus, the process Ijlj — > x\ v m will 
only generate events with 2 leptons + If. At this point, it is necessary to specialize our discussion to a 
single dominant coupling constant hypothesis, assuming Xijk ^ not necessarily identical to X m jj. One 
may distinguish the following four distinct cases. For an LSP x.11 namely, to-^o < mj, (Case 1), only 
the direct RPV three-body decays, x? - * vrfjlk, X? ~* ^djh, are allowed. The branching ratios are then 
determined on the basis of simple combinatoric arguments. For a dominant coupling constant, say, X m ii, 
there are four final states: i^iZ~e + , vil^e - , v m e~e + and D m e + e~ . Accordingly, the branching ratios of 
Xi into two charged leptons will depend on the type (flavor, charge) of the final state: The branching 
ratios equal \ for the flavor diagonal e + e~ or flavor non diagonal Z ± e =F channels, j for the fixed charges 
and flavors l + e~ or l~e + channels and 1 for the lepton-antilepton pairs of unspecified flavors. For a 
dominant coupling constant Ayfe 7^ X m ii, an analogous result is obtained. For m^o > mj, mp (Case 2), 
the branching ratio for Xi decay is, 

3F(xi -> Ijlj) + 3r(xi -> ViVi) 6 
where we have used the fact that in the present case, assuming a dominant coupling constant Ayfc, 
B{lj — > Vilk) = B(pi — > Ijlk) = 1. The factors 3 in the denominator account for the number of families. 
For the intermediate case, mo > m^o > mj (Case 3), there occur contributions from 2-body RPC decays 
and 3-body RPV decays, such that: 

9jlh) = r(x;->^) + r(x^^(^PA) 1 
Kl J r'(x?^p« + 3r(x?^^) 3' y > 
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Figure 9: Branching ratios for the chargino \\ decays as a function of p. The results in the four windows are obtained with 
the following choices for the parameters, [(M2(GeV), mo(GeV), tan/3, \ijk), Tn,[> L (GeV),mj (GeV)]: (a) [(80,20,2,0.05), 
53.19,81.66], (b) [(80,20,50,0.05), 34.20,86.97], (c) [(80,20,2,0.1), 53.19,81.66], (d) [(200,100,2,0.05), 195.6,205.2]. The 
final states are labeled by the letters, A,B,C,D,E, which have the same meaning as in Table M. 
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Figure 10: Branching ratios of the sneutrino decays as a function of fi. The results in the four windows are obtained with 
the following choices for the parameters, [(M2(GeV), mo(GeV), tan/3, \ijk), mp L (GeV),mr (GeV)]: (a) [(80,20,2,0.05), 
53.19,81.66], (b) [(80,20,50,0.05), 34.20,86.97], (c) [(80,20,2,0.1), 53.19,81.66], (d) [(200,100,2,0.05), 195.6,205.2]. The 
final states are labeled by the letters, A,B,C,D,E,F, which have the same meaning as in TablelA. 
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Figure 11: Branching ratios of the slepton decays as a function of fi. The results in the four windows are obtained with 
the following choices for the parameters, [(M2(GeV), mo(GeV), tan/3, \ijk), mp L (GeV),mr (GeV)]: (a) [(80,20,2,0.05), 
53.19,81.66], (b) [(80,20,50,0.05), 34.20,86.97], (c) [(80,20,2,0.1), 53.19,81.66], (d) [(200,100,2,0.05), 195.6,205.2]. The 
final states are labeled by the letters, A,B,C,D,E,F, which have the same meaning as in Table li. 
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Fi gure 12: Mass spectrum of the supersymmetric particles (a), in GeV, and branching ratios for the decays of the chargino 
(b), sneutrino (c) and slepton (d), as a function of fi. The results are obtained for mo = WOGeV , using equation M). 
The final states in figures (b),(c),(d) are labeled by the letters, A,B,... which have the same meaning as in tables MI§Ji 
respectively. 
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where the prime on T' is a reminder to indicate that the decay width includes only the contribution from 
a virtual sneutrino exchange. The approximate equality in eq.(Q) derives from the fact that T'(xi ~ * 
vJjlk) « T(Xi - * ljh)< based on the expectation that an RPV 3-body decay should be much smaller 
than an RPC 2-body decay. An analogous argument to that of case 3 holds for the other intermediate 
case, mj > m^o > mo (Case 4). For the cases 2, 3 and 4, the multiplicity factors are the same as for 

the case 1. The Xi process may occur at the end stage in the decays of Xi ■> I an( i to be discussed 
below. The associated Xi decay multiplicity factors for the two leptons final states will then take the same 
values as quoted above for the various selection criteria. In quoting numerical results below, we shall, for 
convenience, assume the case of unspecified lepton flavor and charge and thus will set the multiplicity 
factors to unity. 

4.2.2 Lowest mass Chargino 

The results in Fig. ^ for the high tan/3 case show a high degree of symmetry with respect to fi <-> —fi, 
which arises from the symmetry in the inos mass spectrum (Figure ||(a)(b)). As can be seen from Figure 
||(a), a dominant mode for the chargino at high values of \/j,\ is the cascade decay, x~ ~ * X°l~v, since 
this occurs via the two-body decay, x~ ~ * (C(7)). Indeed, for these high values of /i, one has 
m§ < vriy- • This two-body decay competes with the other two-body decay E, x~ — > x°W~ , when the 
latter is kinematically allowed, as is the case for // < —200GeV in Figure ||(d). The difference between the 
values of the branching ratios C(7) and E is explained by the relative phase space factors of the associated 
rates. The RPV direct decays (A and B) are three-body decays with small coupling constant and are 
thus suppressed. In the case, < m^o < m^±, the only open channel for the sneutrino is, v — > 11, so 
that the dominant mode for the chargino decay is the RPV decay B(4) (high values of |juj in Figure ||(b)). 
Even for mg ~ m^o, the channel B(4) is competitive due to a small phase space (p ~ — WOGeV in Figure 
||(a)). In this case, for Ay& = 0.1, the direct RPV decay B(4) can become dominant (moderate negative 
values of fi in Figure ^(c)). For small values of the difference between the two dominant leptonic 
(C) and hadronic (D) cascade decays is due to the flavor and color factors. We note also that in a small 
interval of /i near /i = 0, m-yO > m~± (see Figure 0), and consequently the only open channels are the 
direct RPV decays (Figure ||(a)(c)(d)). In this region, the direct RPV decay A(l) is negligible because 
the branching ratio depends on U\\ which is small ]5lj |. In conclusion, the highest branching ratios are 
associated with the cascade decays, C,D and E, except for the case in which the sneutrino is the LSP, 
where they are associated with the RPV decays B. The range of \i for which the chargino x[ is the LSP 
is excluded by the experimental constraints on the inos masses (see Figure ||) . 

4.2.3 Sneutrino 

We turn now our attention to the sneutrino decays. For high values of fi, the cascade decay D has the 
highest probability (Figure |l(](a)(d)) since the decay into chargino is either kinematically forbidden or 
suppressed by a small phase space. As for the chargino study, the RPV direct decay A is of course small 
except when the competitive channel is reduced by a small phase space factor (fi « — lOOGeV in Figure 
|lC](a)). In such a case, the RPV direct decay A may be important for values of near 0.1 (negative [i 
in Figure [l0](c)). When the sneutrino is the LSP, the RPV direct decay has a branching ratio equal to 
unity (Figure |l0|(b)). For small the decays E(6) and F(8) through charginos dominate the decay D 
through neutralinos. The reason is that for [i = 0, Xi is a pure higgsino, whose couplings are weak. In 
the so called higgsino limit, \i — > pi] , the decays B and C(3) are small since they occur through the x~ 
RPV direct decays. However, they have the highest probability if, m^o > m~± (Figure pX|(a) (c) (d)) ■ The 
relation between the leptonic (E) and hadronic (F) cascade decays can be explained in the light of the 
study on the chargino. We conclude that the cascade decays, B, D, E and F, are always the dominant 
modes, except when the sneutrino is the LSP. 

4.2.4 Slepton 

Finally, we concentrate on slepton decays. For high values of fi, the cascade decays via charginos are 
reduced because of a small phase space (\fi\ w 4006^ in Figure pd|(b) or for \i < in Figures pT|(c)(d)) or 
even closed (for fi < in Figure |ll|(a)). In these cases, the decay D via neutralinos dominates. Elsewhere, 
the decays via charginos have higher branching ratios (for /j, > in Figures 0(a) (d)) since larger coupling 
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constants are involved. In the higgsino limit, the slepton cascade decay D via Xi is suppressed for the 
same reason as in the sneutrino study. The decay via \i is then dominating. The interpretation of the 
difference between the decays, B, C, E and F, via charginos is based on the specific behaviours of the 
chargino branching ratios which have already been described above. We see in Figure [ll](c), that for 
Kjk =0.1, the RPV direct decay A is still very reduced. This is due to the important phase space for the 
slepton decay into neutralino. Lastly, a new phenomenon appears for the slepton case. In the higgsino 



limit (/i — > 0) at large tan/3, the matrix element Un — > which forces the vertex I v (see cq.B.3 in 
Appendix |b]) and the branching ratios for the cascade decays through the chargino to vanish . This 
is the explanation of the fact that for /i ~ 0, one observes a peak of the direct RPV decay branching 
fraction (Figure O(b)). Similar peaks are also observed at shifted n < for the low tan/3 cases (Figures 
|ll](a)(c)(d)). However this behaviour appears for ranges of the parameters which are forbidden by the 
bounds on the inos masses (Figure ||). The conclusion is that the cascade decays have always the highest 
probability for the reason that the L-chirality slepton cannot be the LSP in generic supergravity models. 

4.2.5 Discussion 

In summary, we have learned that the general behaviour of branching ratios is mainly determinated by 
the phase space and thus by the ordering of the supersymmetric particles masses. We have explored all 
the characteristic cases, m& > m~- > m^o m~- > > to^o and m-- > to^o > to^. For high values of 
mo lying above M2, the sleptons would have masses greater than the inos masses. We have not analysed 
this case since one has then the same situation in the mass ordering as for the case of small values of 
(except for large enough values of mo where the on-shell W production can take place in I and v decays 
G). In this situation, as we have explained above, the charginos principally decay into neutralinos, while 
the sleptons and sneutrinos decay into charginos. The main conclusion is that the cascade decays are 
the dominant modes except if the sneutrino is the LSP. In this case, the RPV decay, Xi ~ > h&i — > hljh, 
is dominant for the chargino decays, and the only open channel for the sneutrino is of course the direct 
RPV decay. Besides, for values of Xijk higher than 0.05, the RPV direct decay branching ratios can reach 
significant levels for the case where the cascade decays are suppressed due to small phase space factors. 

The excitation of the second neutralino X2 deserves some attention since this may have in certain 
regions of the parameter space comparable, if not larger, production rates than the excitation of Xi- 
Assuming that the direct RPV widths are small enough so that the decay chain is initiated by the RPC 
contributions, then the desintegration mode, x° — * (Xi + (Xi + vv), will also yield 21+i? and 41+^ 

final states, respectively, and the other desintegration modes, x\ ~ ¥ (xt X\ (pv, viz), {I P), 

will yield 21+^ and 41+J^ final states according to decay schemes similar to those given in Tables 
In our supergravity models, the X2 decay into Xi should be suppressed by a small phase space (Fig.[3|). 
To determine which of the decay modes, X2 ~ ¥ Xi? I 01 v-> leads to the dominant signal would require a 
detailed comparison of branching ratios at the initial as well as the subsequent stages. 

Let us ask in what way would alternate hypotheses on the family dependence affect our conclusions. 
Especially regarding the multiplicities of final states, this is relevant for the cases, m^- > mj, mp, where 
the chargino can cascade decay to on-shell sle pton s or sneutrinos (A(2) and B(4) in Table [I]). As we have 



emphasized in the last paragraph of Section 2.2, the chargino decays have a multiplicity of 2 for three 
degenerate families of sleptons. For the case of two degenerate families, labeled by the indices, to, n, 
assuming a dominant RPV coupling constant Xijk, the multiplicity equals 2 for (m, n) = (i,j), since the 
two sleptons from families i and j can be produced on-shell, and equals 1 for m — k or n — k. For 
the physically interesting case of a single low mass family, labeled by the index, to, one finds that the 
multiplicity equals 1 for to ^ k and otherwise. The conclusion is that the RPV contributions A and B 
(in Table to the chargino branching ratios increase as the number of slepton families, which are lower 
in mass than the chargino, becomes higher. This effect, which is quite small, would affect the branching 
ratios in parameters regions for which the RPV contributions A and B are not weak, that is for fi < in 
Figures p|(a)(b)(c). 

In Figure we present results for the branching fractions for fixed too in the infrared fixed point 
model with electroweak symmetry breaking. In this constrained version, where mi/ 2 varies with the 
dependence on fi is rather similar to that of the non minimal model where we worked instead with fixed 
mo and m^. However, as we see from the mass spectrum, here the LSP is the neutralino Xi f° r a U the 
physical ranges of the parameters. Due to the large mass difference between the x? LSP and the NLSP 
(next to LSP), the cascade decays are the only dominant modes and the branching ratios for the RPV 
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direct decays are very weak. 

Let us add a few qualitative remarks on the predictions of gauge mediated supersymmetry breaking 
models. In order for the production rates in the minimal model [5^ ] to have the same order of magnitude 
as those obtained in the supergravity model of section |I| one needs a parameter A = -p- ~ l0 4 GeV, using 
familiar notations for the supersymmetry breaking scale (y/~F) and messenger scale (M). Concerning the 
signals, by comparing the mean free paths for Xi (favourite candidate for LSP) in both models, one finds 
that the decay channel to the gravitino, xi ~~ * 7 G, becomes competitive with the RPV decay channel, 
y9 -> vll for V<F> < iSll 

Let us also comment briefly on some of the experimental issues. A given final state can possibly arise si- 
multaneously from several of the single production processes. The important 41+^? signal which occurs for 
X , I , v productions is one such example where one may be forced to add all three types of cross sections 
in comparing with some given experimental data sample. Similarly, for most signals, one must typically 
add the contributions from the two charge conjugate partner processes. Concerning the competition with 
the standard model background, one expects that the most important contributions to the final states, 
21+^ and 41, will arise from the reactions, ijlj -► W + l~v, W~l + v, W+W~ , Z°l + l~ , Z°Z°, Z° 7 . In 
spite of the large standard model rates of order one picobarn at s/s = 500GeV , one should be able to 
distinguish the single production signals by exploiting their specific non diagonal flavor character (final 
state B in Tableland A in Table ||). The other multileptons final states, generated by the cascade decays, 
41+i^, 4l+Z°, 3l+Z°+W ± +I$,... have a standard model background which is negligible. The potentially 
large two photons background processes, induced by 77 photons pairs radiated by the initial leptons, 
can be significantly reduced by imposing suitable cuts on the leptons transverse momenta. Finally, we 
note that the selection by the RPV single production of identical helicities for the initial state, can 
be exploited to discriminate against the minimal supersymetric standard model and also the standard 
model, for which the identical helicities configuration only appears with the t-channel Z-boson exchange. 



5 Dynamical distributions 

The distributions of rates with respect to kinematical variables associated with the final states offer 
helpful means to characterize the underlying production processes. As an indicative study we shall 
present here some characteristic dynamical distributions obtained for the production reactions, l^jlj — * 
Xi P, Xi v i X?£j X% v i X^Pi from a Monte Carlo events simulation for which we have used the event generator 
SUSYGEN Q. We concentrate on the final state signals of 21+$, 41 and 41+$. Note that for high values 
of /1, the final state 41+i? is the dominant mode for the chargino, slepton and sneutrino decays. This 
signal also receives contributions from the reactions, l^lj — > DZ (DZ°), 2/7(1/7), / ± M /=F , which however 
are not included in the simulation. The standard model background is expected to be small for the 
41+_^ signal. The main background from the minimal supersymmetric standard model interactions arises 
from the neutralino RPC pair production, — > XiXi- Following the analysis in |33|, we consider an 
incident energy of i/s — 350GeV and use a non minimal supergravity model for which we choose the set 
of parameters, M2 = 250GeV, [i — 4Q0GeV, mo — 70GeV, tan/3 = 2, which yields the spectrum, m^o = 
118.5GeF,m^o = 221.4GeV,m-± = 219.1GeV,m«> r = 22bGeV,mr = 233GeV,m r = UIGeV. The 

integrated rates (ignoring acceptance cuts) are, for \ m jj = 0.05 [m — 2], <t(xiH~) = 30.9/6, a(xi>^^) = 
4.8/6, cr{x2^ij) = 12.1/6 and cr(xiXi) = 238.9/6. We consider the following five dynamical variables for 
all types of final states: Invariant missing energy, E m = ^ iei/ Ei where the sum is over the neutrinos, as 
appropriate to a broken R parity situation; Average per event of the fi lepton transverse momentum, 

P^/i*) = ^ ; where is the number of muons; Angle between the momenta of same electric 

charge sign (SS) muons pairs; Average per event of the summed transverse momenta for leptons pairs 

of same sign (SS) or opposite sign (OS), Pf S ' OS (ll) = Yluj) — n'^^ 3 ' ' wnere N is the number of 
configurations and I — e, /1. 

We have generated the inos single production and the Xi P a i r production in separate samples of 5000 
events each. Our choice of using equal number of events for both reactions has been made on the basis 
of the following three somewhat qualitative considerations, none of which is compelling. First, the single 
production reactions occur in company of their charge conjugate partners, which multiplies rates by a 
factor 2. Second, the other v and I single production reactions, which have not been included, would be 
expected to add contributions of similar size to the leptonic distributions. Third, assuming for the RPV 
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Figure 13: Distributions of missing energy, muon transverse momentum, same sign muon pair angle, summed transverse 
momentum for opposite sign (OS) and same sign (SS) leptons pairs (electrons and muons) for the single production 
processes, ij/J — + xt > Xi^mj Xi^Vi X-j^M' X-jiV (solid line), and the pair production process, JjJJ — * X?X? (dashed line), 
at a center of mass energy of 350G?eV. 77ie parameters values are, M2 = 250GeV, mo = 70GeV, fi = AOOGeV, tan/3 = 
2, A211 = 0.05. Events samples, consisting of 5000 events each, are generated for the inos single production and neutralino 
pair production, respectively. 



25 



coupling constant the alternative bound, \ m jj 100GeV < 0.05, there would result a relative enhancement 

for single production over pair production by a factor, ( ioogIv 7 ^ ne & bove three points motivate our 
rough guess that the number of events chosen for the five single production processes (together with their 
charge conjugate partners) should be of comparable size to that of the Xi P arr production process. 

The results are shown in Figure [13[ The single production reactions present certain clear characteristic 
features: Concentration of missing energy E m at low energies; pronounced peaks in the muon transverse 
momentum P t (/n ) and in the angular distribution for the same sign muons pairs; and a double peak in 
the transverse momentum distribution for the opposite sign leptons pairs, Pf s (11). The large transverse 
momentum components present in the single production distributions in Pf s (ll) and in P t (/i ± ) are 
explained by the fact that one of the two leptons (namely, V^) is created at the production stage. Similarly, 
the existence of a strong angular correlation between same sign muons pairs is interpreted naturally by 
the momentum conservation balance between the lepton 1^ produced in the initial stage and the other 
lepton produced at the decay stage. Although there are certain distinguishing properties between the 
single and pair production processes, the discrimination between the two may depend crucially on the 
relative sizes of the associated event samples. Of course, the best possible situation would be for an 
energetically forbidden neutralino pair production. 

Finally, we comment on the effect of eventually excluding the single production component. In that 
case, most of the signals for single production would become less diluted in comparison with neutralino 
pair production, the large missing energy signal would be removed while the large OS lepton pair Pf s (11) 
signal would become amplified. 

6 Conclusions 

We have analysed the full set of 2 — * 2 single production processes at leptonic colliders induced by the RPV 
interactions LLE C , within a supergravity model. Although our approximate study has obvious limitations 
(factorisation and narrow resonance approximation, neglect of the spin correlations and omission of 
acceptance cuts), it uncovers the general trends of all the 5 single production reactions. Over the whole 
parameter space, for an RPV coupling constant X m jj of order 0.05, the integrated rates are of comparable 
order of magnitudes although X\ an d v are typically larger by factors of 2 compared to I production and 
by factors of 5 compared to Xi production. The detectability for each single production separately is 
modest at LEPII but comfortable at NLC, corresponding to a few events and a few thousands of events 
per year, respectively. A wide region of the parameter space can be probed at \ m jj > 0.05. In spite of 
the rich variety of final states, the dominant signals arise from the single or double RPC induced cascade 
decays to the LSP, which are also favored by phase space arguments. For the minimal supergravity model, 
assuming electroweak symmetry breaking, large mass differences occur between the scalar superpartners 
and the Xi LSP, leading to dominant cascade decays modes with weak competitivity from the RPV direct 
decays. The signals, 41+^,61,61+^?, arising from cascade decays are free from standard model background 
which make them quite interesting signatures for the discovery of supersymmetry and R parity violation. 
Even the 41 signal arising from direct RPV decays could be observable due to a characteristic non diagonal 
flavor configuration. For center of mass energies well above all the thresholds, the 41+ J^? signal receives 
contributions from all five single production processes and hence should be strongly amplified. We have 
presented some dynamical distributions for the final states which could characterise the single production 
reactions. 
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A Formulas for spin summed amplitudes 

We discuss the five 2^2 body single production processes given by eq. ([!]). The formulas for the 
probability amplitudes are: 
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In deriving the results for the inos production amplitudes, we have systematically neglected their hig- 
gsino components. The parameters in the Z°ff and Z°ff vertices denoted as, a#(/) = a(fii) and 
o-nif) = a (fii), are defined by, a(fii) — o-ifn) — 2T^ (/) — 2Qxw, with H = [L, R] and xw = sin 2 9w- 
Throughout this work, our notations follow closely the Haber-Kane conventions [[l8| . 

The unpolarized cross sections in the center of mass frame are given by the familiar formula, da j d cos 6 = 
p/(1287rfcs) J2 P oi \ M ?> where the 

sums over polarizations for the probability amplitudes squared are given 

by: 
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where Re stands for the real part, R s (vi) = s — m 2 . + im^Tp, Rt{vi) = t — m 2 . and R u (vi) = 
u — m 2 ,., [s = (k + k') 2 , t = (k — p) 2 , u = (k — p') 2 ], with similar definitions applying for the propagator 

factors R s ,t,u(h,h)- 

B Formulas for partial decay widths 

The formulas for the various two-body decay widths are quoted below. 
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We use the notations: 



m~±, M 



m^o and fc 



\i(M 2 ,m 2 ,m%)/2M with X(a,b,c 



O u = 0« m = N l2 U ml 



a 2 + b 2 



-ab- 



M± = 
bc + ac. The notations, 



X"/, Qf , stand for the third component of the SU(2)l group and the electric charge of the fermion /. We 
have omitted the higgsino components of the inos. We shall use the simplified formulas for the RPC three- 
body decays, x~ — ► X? + ^ <7<L obtained by neglecting the three-momenta in the W and / propagators, 
as quoted in |54}] . We have set in these formulas, the flavor and color parameters to, Nf = 2, N c = 3 for 
quarks and Nf = 3, N c — 1 for leptons. The formulas for the spin summed amplitudes of the RPV decays 
Xa ~ ¥ ViVjlj~, Xa ~ > ^t^JK i associated to the coupling constants Xijk, were first derived in |l5j (see 
the appendix). The integrated decay rates are given by familiar formulas p7j | involving twofold integrals 
over the final state three-body phase space. If we neglect the final particles masses, an analytic formula 
can be derived for the integral giving the contributions to the charginos partial rates associated with 
the gauginos components only (neglecting the higgsino components contribution). For completeness, we 
display the final results: 
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where Sp(x) — Polylog(x) — Li 2 (x) is the Spence or Polylog function. We use the notations fi a = 
ml JM 2 , [a = X al = U al for the decay Xa ~> l t l 7 l i '> and Ma = m? /M?_, [a = i, j], X a i = Ki 

Xa J l a Xa 

for the decay x~ ~~ * ViVjlT 
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